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ABSTRACT 


In this paper the review of mathematical modeling of two phase blood flow in capillary is done with disease Diabetes .There are 
two phase of blood flow, one is RBC and other one is plasma. By the Fahreaus-Lindqvist influence the blood flow in two sepa¬ 
rated layers as transient from side to side capillaries, in which one layer is Newtonian and other is consider as Non-Newtonian 
layer. Further Non-Newtonian power law model is applied to bio fluid mechanical system. We have collected an observation data 
for patient having Diabetes. On the whole arrangement is in tensorial form and key method adapted is investigative as well as 
arithmetical. Later the graphical representation is also given for the data used here. 
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INTRODUCTION 

Composition of Kidney: Kidneys are bean shaped organs. 
Kidney has an outer fibrous renal capsule and supported by 
adipose tissue. There are two main parts, inner medulla and 
cortex. External cortex is reddish brown, where fluid is clean 
from blood. Inner medulla is paler and prepared by coni¬ 
cal formed sections call renal pyramids. [1][2] The medial 
edge of the kidney is called the hilus and is the part where 
the renal blood vessels depart and come in the kidney [3][4] 
[5] [6]. The main purpose of kidneys is giving out the blood 
and takes away waste and surplus water through the urine. 
Both kidneys are comprised of million filter systems named 
nephrons. Each nephron filters little quantity of blood. It has 
a filter that is glomerulus, and a tubule. Glomerulus allow 
fluid and waste items bypass through it, although, it prevent 
blood cells and big molecules, typically proteins, from pass¬ 
ing. The clean fluid then passes from tubule, which sends out 
essential minerals reverse to the blood stream and eliminates 
waste. [7] [8] 

Blood contribute to kidney: In the stomach, the renal artery 
stem from abdominal aorta lower to the bigger mesenteric 
artery and enlarge crossways in the direction of the kidneys. 


Just earlier to reach the kidney, all renal artery separates keen 
on five segmental artery, which supply blood to the variety 
of regions of kidney. Every segmental artery goes into the 
hilus of the kidney furthermore divides into quite a few in¬ 
terlobar arteries; get in front of the renal columns linking 
the renal pyramids and take blood in the direction of the pe¬ 
ripheral of the kidney. At the link in the midst of the cortex 
and medulla, the interlobar arteries figure out arcuate arter¬ 
ies, which spin to go behind the contours of renal pyramids. 
Arcuate arteries form a number of branches, identified the 
same interlobular arteries, divide at right angles and enlarge 
the entire mode in the course of the renal cortex on the way 
to the outside of the kidney. In latent fully developed kidney 
get 1.2 to 1.3 1 blood per minute. The Fick principle is ap¬ 
propriate to find renal blood flow [9] the renal blood flow 
is considered by dividing by one minus the hematocrit. [10] 

Strain in renal vessel: Blood is a fluid tissue consists of ap¬ 
proximately 55% of fluid plasma and 45% of cells. It consists 
of three major types of cells which are red blood cells, white 
blood cells and platelets. Blood plasma is made of 92% wa¬ 
ter and 8% of ions, Proteins, metabolites. The standard thick¬ 
ness of whole blood for a human is on 1060 kg/M 3 .[ll] The 
force in glomerular capillary has been calculated directly in 
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the rate and has been set up to be significantly lower than the 
predicted on the source of not direct dimension . When the 
mean systolic arterial force is 100 mmhg, then glomerular 
capillary force is about 45 mm hg. The pressure drop cross- 
ways the glomerulas is only 1 to 3 mmhg , but supplementary 
jump down occurs in the efferent arteriole such that pressure 
in the per tubular capillary as regards 8 mm hg. Pressure in 
renal vein is regarding 4 mm hg. [9] 

Structures of renal capillary: Renal vascular model is ex¬ 
traordinary in this blood flows through two capillary beds, 
high pressure (glomerular) and low pressure (per tubular), 
associated in sequence. Blood enters the kidney by means of 
the renal artery and, following a chain of divisions, arrives 
at the glomerulus. Glomerular capillaries have to first go by 
afferent arteriole, where blood passes all the way through a 
second arteriole, the efferent arteriole. After that blood flows 
through the per tubular capillaries, which comprise the vasa 
recta that enlarge into the renal medulla, after this it accumu¬ 
lates in gradually larger venules and veins, and then goes out 
of the kidney through the renal vein. [6] [12] [13] 


MATERIAL AND METHODS 

Disease (Diabetes): Diabetes is a disease in which the body 
does not correctly absorb food for make use of energy. The 
majority of the food we eat is converted into glucose, for our 
body to utilize for energy. The pancreas makes a hormone 
named insulin to assist glucose for entering the cells of body. 
Body having diabetes will some time doesn’t make sufficient 
insulin or can’t take up its own insulin as how it should. This 
effect sugars to put up in blood. So diabetes is also referred 
as sugar. [14] [15] 

Explanation of the problem: Capillaries are very thin 
as well as distant from the heart, so in this situation how 
the blood flow is possible in these vessels. This is done by 
Fahreaus-Lindqvist effect. As per this effect the blood flows 
in two alienated layers while transitory through capillaries. 
The plasma layer contains more or less no blood cells. The 
second layer is blood cells which hang in plasma on the 
alignment of the capillary. In this procedure the useful blood 
viscosity depends upon radius of the capillary. So the effec¬ 
tive viscosity decreases, with the radius and thus the blood 
flow becomes probable [6. 

Model: Blood viscosity means thickness and tackiness of 
blood. Usually mature blood viscosity is 40/100, unit mill 
poise.[16] Blood is a dynamic organ in so far since it act as a 
non-Newtonian fluid, in the sense its viscosity changes with 
role of shear rate. Consider shear rate like velocity, when 
blood travels fast as in peak-systole, it is physically thin¬ 
ner; when it travels slowly during end-diastole, it is thicker 
and stickier, because red cells cumulate. This method is 
called as the shear-thinning, non-Newtonian character of 


whole blood. [17] [18] [19] [20]. Here in the current paper we 
choose comprehensive 3- dimension orthogonal curvilinear 
co-ordinate arrangement, which is arranged as E3 known as 
three-dimensional Euclidean space. Here are some quantities 
related to moving blood in cylindrical vessels: blood velocity 
V k = V k (x\t) ? k=l,2,3 blood pressure P = p(x i ,t) and density 
p=p(x i ,t) where x 1 is co-ordinates of any random point in 
space and i-1,2,3. 

Hematocrit: It is the relative amount of the quantity of red 
cells to the quantity of whole blood. Normal array for hema¬ 
tocrit is dissimilar among the sexes and is roughly 44% to 
51% for men and 36% to 47% for women. The hematocrit 
(articulated as percentage) is usually as regards three times 
the hemoglobin (grams per deciliter). It is denoted by H.[6] 
[21] [22] [23] [24] 

Equation of Continuity: As there is no source or sink in 
the whole circuit of the human blood circulatory system, the 
heart behaves simply like a pumping station that is why the 
law of conservation of mass is applied to hemodynamic [25]. 
In view of the fact that, entire blood flow circuit of the kidney 
is known as a Renal Circulatory System. Consequently renal 
circulatory system is a subordinate system of human circula¬ 
tory system. Blood come into kidney by arteries and out by 
veins and in a kidney no starting place or is submerged. 

Mass of enter the blood = mass of outer the blood. Hence 
law of conservation of mass is applied for renal circulatory 
system. The blood flow exaggerated in presence of blood 
cells. This outcome is directly relative to volume engaged 
by blood cells. Let X is the volume portion enclosed by the 
blood cells in unit volume. And X can be taken as H/100. 
Therefore volume portion of plasma is (1-X), 

If mass ratio of blood cells to plasma is r, then we have 

r = —^— Everywhere p c and P p are densities of blood 
( 1 - X )P P 

cells and blood plasma. Actually above mass ratio is not 
steady; even then it may be hypothetical to be constant in 
current situation. [26]. 

The two phase of blood, i.e., blood cells and plasma go with 
an ordinary velocity. In this paper we have used the model 
given by Campbell and Pitcher. As per this model we consid¬ 
er the two phases of blood separately [27]. As per the princi¬ 
ple of conservation of mass, the equations of continuity for 
the two phases are as follows [28]. 

U(Xp c V‘),i = 0 

X) p n 

^+ ((1-X)p p vy = 0 

Here v is common velocity of two phases and (X/^VAi as co¬ 
variant derivative of (X/^Vyi with respect to X'. As a result 

■ 



9{Xp, 
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((l-X)p V with respect to X 1 If we denote uniform density p m 
by: 


1 + r _ r + 1 
Pm Pc P$ 


Equations can be collective together as follows, 


dt 


+ (y 0 m c),i = 0 


( 1 ) 


As we know that blood is incompressible liquid hence p m 
will be a content measure. Hence the equation of continuity for 
blood flow takes the following form: 

V,| = 0 

av 1 | v'd'/g i Q 

ax J,gx' 

Equation of Motion: As per this principle, total momentum 
of every fluid system is preserved in lack of outside force. So 
the law of conservation of momentum is appropriate to renal 
circulatory system. In addition, rate of change of momen¬ 
tum of a fluid particle with respect to time equals to external 
force exerted on it, which is also known as Newton’s second 
law of motion.Therefore, rate of change of momentum as 
same as sum of regarding two mentioned forces, which could 
be symbolically presented as follows, 
dp dp 

= -P+F where, = rate of change of momentum, 
P=Intemal pressures=viscous force 

The hydro dynamical pressure p among phases of blood could 
be hypothetical to be consistent since phases i.e. blood cells 
and plasma is for all time in the symmetry state in blood [29]. 
Consider viscosity coefficient of blood cells to be Tl c , apply¬ 
ing the principle of conservation of momentum, we conclude 
the eq. of motion for phase of blood cells: 

Xp c ^T + (xp c v i )v j j = —Xp,j g L > + XT] c (gi k v4) ;j (2) 

Take viscosity coefficient for plasma as n^. Equation of mo¬ 
tion of plasma is: 

(1 “ x)p P^ + ( (1_x)p P vi ) v 'i = -ft- K ^e li + ( 1 - x )ip(g il ‘<)i' (3) 

On adding up eq. (2) and (3), put in relation (1), equation of 
flow of blood with both phases is: 

P m ^r + (pm vi N = g Li + v(g ik vi)^j 

Where n m = Xn c + (l - X) ii p i s viscosity coefficient for blood 
as a combination of two phases. 

Special constitutive equations for blood: Normally blood 
is non-homogeneous combination of plasma and blood cells. 


Although for practical reasons it can be considered to be ho¬ 
mogeneous two-phase mixture of plasma and blood cells. 
The constitutive equations planned for whole blood mixture 
are as follows: 

(i) Newtonian equation: 


x=r|e, 

Where q is the viscosity coefficient 

This hold fine in the broad blood vessels where there is low 
hematocrit [30]. 

(ii) The non-Newtonian power law equation: 

e° 

x=q 

This is conformable for strain rate between 5 and 200; 
0.68<n< 0.80 [10] 

The non -Newtonian Herschel-Bulkley equation [10] 

x=1] eV°(T> TO) 

e=0(r < Tq). It holds superior when blood shows yield stress 
Tq. we notice that the yield stress arise because blood cells 
form aggregates in the form of rouleaux at low strain rate. 

If t < tq ? no blood flow takes place. It is found that yield 
stress is given by the following formula 

| A(H - H m ) 

Tq 100 

1 

Where, A= (.008 ± Q.002dyite/cm 2 >; H m is the hematocrit be¬ 
low which there is no yield stress. 

Boundary Conditions: (i) The velocity of blood flow on the 
axis of capillaries at r=0 will be utmost and finite, say V Q 
= maximum velocity. (ii)The velocity of blood flow on the 
fence blood vessel at r=R, where, R is the radius of capillary, 
will be zero. This state is well known as no-slip condition. 

Mathematical Modeling: Consider the two layer blood 
flow in that one is Newtonian while other is Non- Newto¬ 
nian power law flow. The layer which is close with wall of 
the vessels can be taken as Newtonian; reason is layer con¬ 
tains plasma only. The second core layer can be considering 
as non-Newtonian power law, reason here the ratio of blood 
cells is too high in comparison to plasma. 

Equation of continuity for power law flow will be: 


y/s[y/s^) i (4) 

Again the equation of the motion is extended as: 


Where T^ is taken from constitutive equation of power law 
flow 
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Pm % + An v ' v ’j 1=T ’j' j 


(5) 


= -pgV + rj m (Jj) n =-pgV + T V And Pm~ x Pc + ( 1 ~ x )Pp 

n m =Xri c +(\-X)ri p 

Because blood vessels are cylindrical, the above equations 
should transform in cylindrical co-ordinates system. Now we 
have to transform equation (4) and (5) in cylindrical form. 

Cylindrical Co-ordinates, X 1 =r, X 2 =0 , X 3 =z 

Matrix for metric tensor as cylindrical co-ordinates: 

"1 0 0' 

0 r 2 0 

0 0 1 


T 'i' ij=T ' ij)+ • U T ' ij 

J y/SdX J [j k J 

Keeping in view the above fact, the governing tensorial 
equation can be distorted into cylindrical form : 


Equation of continuity - 


dv 

dz 


= 0 


g ij 


Equation of motion 
R-component 


Matrix of conjugate, metric tensor: 

n o o 

aT 


- 3 e=o 

dr 

6 - component 
0 = 0 


g~ 


0 l/r 2 0 
0 0 1 
Christoffel’s symbols of 2 nd kind as: 

f 1 ] f 2 1 f 2 l i 

(2 2j [2 if }l 2( r 

Relation among contra variant physical components of ve¬ 
locity of flow of blood is: 

Vgll vl= v r = 

•v/§ 22 v2= v 0 


Z-component 


n= 

0z r 5r 


0V Z 

Sr 


n 


(9) 


Considering flow of blood as axially symmetric in arteries 
i.e. V =0 and V 

u r 

V z and p do not depend upon 0. Also the blood flow steadily, 
i.e. 


> v r = v* 

2 

Vq= rv z 


Sp_Sv r _ v fl _Sv 7 _ 0 
St St St St 


( 10 ) 


And V 833 V 3 = v z ^v z = v 3 
Also physical components of 


SOLUTION: 

Integrating equation (6) we get, v z =v(r) because v does not 
depend upon 0. (11) 


P :g 1J IS 

?J 


/g H ¥ 


>J 


The matrix of the physical components of shearing stress- 
tensor: 


Integrating equation of motion (7) yields: 
P=p (z) since p does not depend upon 0 


( 12 ) 


T ij - r hn = ! hn (g lk P t + g jk P h f is 

0 0 

0 0 

t] m (dv/ dr) n 0 


k ' s y ,k 
rj m (dv/dr) n 

0 

0 


Now, with the help of equation (11) and (12) the equations of 
motion (9) convert: 


q- _ dp + Um d 
dz r dr 


dv 
dr , 


(13) 


The pressure gradient (dp/dz) = p of blood flow in the arter¬ 
ies remote the heart which is supposed to be constant and 
hence the equation (13) converts: 


The covarint derivative of T 13 


d_ 

dr 


dv 
dr , 


= ._PL 


(14) 


*7i 


m 
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On integrating equation (14), 


dv 

dr 


n 9 

= 

^■ 7 lm 


(15) 


Because velocity of the blood flow on the axis of cylindrical 
arteries is maximum and constant. So apply the boundary 
condition at r=0, v=V 0 (constant), equation (15) converts: 


r 





dv 

dr 


2 

pr^ 

277 m 


1 In 


pr 

2 ; 7m 


(16) 


On integrating equation (15), 


v = - 


nl/n 


2 77i 


m 


^.1/n+l 

n+l/n 


+B 


(17) 


To determine the arbitrary constant B, apply the non-slip 
condition on the inner wall of the arteries at r=R, V=0, where 
R= radius of vessel, on equation (17), 


B= 


1/n ' 

nR ,+n 


2/7, 


m 


(« + !) 


Hence the equation (17) converts: 


V= 


p 

n 

277 m 

n +1 


1 

Rn + 1 



Which conclude the velocity of the blood flow in the artery 
remote from heart. Now the formula for velocity of blood 
flows can be obtained by replacing with in Newtonian 
model: 



R-S < r < R 


Where 5 the radius of core layer. The velocity of core layer is 
obtained as the equation (15) of power law model: 


m 


1 

n n 


2r| 


m; 


+ 


4r( r 


n +1 
2 


R ” —r " 


r2-(R-8) 


1 


P ' 

n n 

, 2T lm, 

n +1 

V 


1+1 1 +] 
Rn -(R-8)n 1 


:0 < r < R-8 


Where, the 2 nd term is the relative velocity of plasma layer 
with respect to core layer. 


STATISTICAL METHODS 


Observation data: 


S.no. 

H.B. 

Hematocrit 

Blood Pressure(mmhg) 

i 

14.9 

44*7 

130/80=17331.6/10665.6p 

2 

14.8 

44-4 

120/80=15998.4/10665.6p 

3 

15.9 

47-7 

130/80=17331.6/10665.6p 

4 

14.7 

44.1 

120/85=15998.4/11332.2p 


Bio-Physical Interpretation: 

The blood flow (Flow flux) in capillary is 



27trdr 


Average systolic Pressure =125 mmhg 
Average Diastolic pressure = 81.25 mmhg 


Pressure at Capillary = 


D+S+d 


= 61.46 


2 2 

Pressure on Venuies= — of Capillary = —><61.46=41 mmhg 
p. =blood /?ressureoncapillary=61.46xl33.322=8193.97ps 
p^ = Mood pressure on venules =41><133.322=5466.2ps 


Now we have, 

q m = viscosity of mixture =3.5 x 1 0"^ ps 
r|p = viscosity of plasma =1.2x10"^ ps 
R= radius of capillary =0.0965m 
5 = thickness of RBC layer =iR=0.0322 
R-S =0.0643 

averageHB = 15.07 & averageHematocrit=3x15.07=45.21 
Q =flow flux of blood = 0.01833 pa 
Length of capillary = 19000m 
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n m =n c x+(i-x)np ; x=i=o. 452 i 


3.5x 10' 3 =r| c xO.4521 +(l-0.4521)x 1.2xl<)- 3 
r( c =0.0063 

dp 

P ~~dz 

=>pjdz = -/dp =>p(z f -z i ] = [p i -p f ] 


Now 


_ Pj-Pf 

z f z i 

4p 


_ 2727.77 

length of capillary 19000 


=0.1436 


R-5 

Flow flux Q= j 

0 


2r| 


r'im; 


n 

~n + \ 


±+i i +1 

-r 11 


41 R "' (R ' 5r 


1 


p 

n JL_ 


n + 1 

V 


*■ -(R-8)n 


UP 


2jtrdr 


+ J -^f+T 2 ) 

R-5 4r lp ' 


2nrdr 


0.0643 

Q = l 


0.1436 W n 
2x3.5x10+ l n+1 


n+1 n+ 1+1 

0.0965 n 2nrdr-2nr n dr 


°- 14 — , f0.0965 2 -0.0643 2 Wdr 
4xl.2xlO‘ J J 1 ' > 

_n_ _n_ 

0.0965 n+1 -0.0643 n +l 


1 


0.1436 V / n 
n +1 


2x3.5x10" 


27irdr 


0.0965—0J436 f 0 .0965 2 2jtrdr-2jtr 3 dr 
0.0643 4x1.2x1 O' 31 ' 


_ ( Yl \ ii~r 1 

Q=(20.5l)»x| —- |x0.0965 n 2 tc 


n+1 


0.0643 


Jo 


-2k 


' 3n+l 

0.0643 


r n 

+0.97 

r 2 

3/7 + 1 

2 

n 

0 



0.0643 


Q=(20.5l)»l — 


(0.0965)” ;r (0.0643) - 2k | |(0.0643) » 


1 -(0.002) -(20.51)-x 


n+1 


(0.0965)+ -(0.0643)" 


0.0129+0.00125 


Q=(20.51)fe 


+ (20.5l)nx__(0.0643)« 


(0.0965)« +1 0.0129-(0.0965)« +1 0.0129 
0.0017n 


0.00083 -- 


3n+l 


+0.00325 


Q—0.01833— (20.51x0.0643)« 


0.00083- 


0.0017+ 


3n+l 


+0.0033 


0.01833-0.0033=(l.32)n; 


n +1 


0.00083- 


0.0017n 


3n+l 


0.01503=(l.32)nx_^_ 

v ’ n + 1 

n=0.045; Again, 


0.00083- 


0.0017n 

3n+l 


R-5 

Flow flux Q= J 
0 


2r l; 


mj 


n 

~n +T 


41 r2 - (r - 8) 


R" + '-(r)n + 1 

2 


n 

n +1 


2 %iJ 

+ f -f£-fj? 2 -r 2> )2jtrdr 
R-8 4l3 p x > 

J-j x( 0 .043){(0.0965)“-r !! |2radr 


R- -(R-S)n 


i+1 


27irdr 


0.0643 

Q= f 

o 


4x1.2x10'- 


0.0965 2 -0.0643 2 )^27irdr 


(0.043)x(0.0965 232 -0.0643 232 ) [27rrdr 


0.0965f 0 

+ i —^ 

0.0643V 4x1.2x10"-' 


0.0965^ 


x2ra-dr-27ir 2 dr| 


/ \22.2 


r r 2 n 

0.0643 

„25.2 ' 

0.0643 “ 

ro 

O 

O 

1 

0.0965 23 ' 2 2ji 

r 

~2 

-271 

0 

r 

25.2 

0 


0.0048 

P 

2f i m 


(0.0052)2^- 


0.0643 



r r 2i 

J27T 

r 

2 


0.0643 


-(20.5 l)nf 


I+-1 

' n+1 n+1 ^ 


r r 2j 

0.0643 


0.058 

r 2 

0.0965' 

0.0965 n -0.0643 n 

2n 

r 





U+u 



2 

0 

+ p 


2 

0.0643 


+29.92 


0.0585 

- 2 J 4 


9 -.0.0965 

2 Jo.0643 

0.0965 

'0.0643 


0.0048 


-271 


,0.0965 


0.0643 


Q= ^— (2.5027x10' 3 ')+0.014p+0.0087p 

l 2r lm J 

Q=|++) (2.5027 xio _3, )+0.0227p 
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Ignoring the term 0.0227p (negligible value) 

( V 2 ' 2 

Q=0.01833= yE— (2.5027X10- 31 ) 

J 


4^— H^^xlO 2 *) 22 - 2 

2l lmJ [ ’ 


as; p : 


Ap 

19000 

Ap 


380001 ^ 


19.96 


we get the relation Ap = 38.68H + 910.176 and the graph 
is drawn between Blood Pressure and Hematocrit in re¬ 
nal capillary in Non-Newtonian flow, and trend of graph 
shows the relation between Blood Pressure drop and He¬ 
matocrit as linear as y = 5.803x + 2645 .This linear rela¬ 
tion approves the two phase relation = ^n c + (1" 3Q \ 
Where X=H/100 and slope of trend line is.... 5.803. By this 
slope of the trend line we can suggest about blood pressure 
fluctuation in regarding tissue. 
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In this paper we have taken observed data with respect to 
Blood pressure and Hematocrit considering the disease Dia¬ 
betes. On applying the two phase Non- Newtonian power 
law model, we get the relation Ap = 38.68H +910.176. 
Later the graph is interpreted, where the trend line shows a 
linear relationship between pressure drop and hematocrit as 
y= 5.803x+2645, where the slope of the line is 5.803. 

CONCLUSION 

In Bio Physical interpretation, we have taken observed 
data regarding with Blood Pressure and Hematocrit of 
Diabetes patient. Using two phase Non-Newtonian model 
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